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Abstract 

We propose a generalization of the quantum entropy power inequality involving condi- 
tional entropies. For the special case of Gaussian states, we give a proof based on perturba- 
tion theory for symplectic spectra. We discuss some implications for entanglement-assisted 
classical communication over additive bosonic noise channels. 

1 Classical and quantum entropy-power inequalities 



s 

S The entropy power inequality, proposed by Shannon j27j and later established with increasing 

Q-( rigor by Stam [2S] and Blachman [3], has become a fundamental tool in classical informa- 

tion theory. Shannon's original application of the entropy power inequality is a lower bound 
'""' on the capacity of an additive (but potentially non-Gaussian) noise channel [271 Theorem 

,__! 18]. However, the usefulness of the entropy power inequality is especially evident in multi- 

C£} terminal information theory. Among the most well-known applications are the characterization 

of the Gaussian broadcast channel by Bergman |1], the Gaussian two-description problem by 
Ozarow [25] and the quadratic Gaussian CEO problem by Oohama [21]. In these multi-user 
O settings, Fano's inequality by itself is insufficient to characterize different tradeoffs. A more 

^^ recent application proposed by Liu and Viswanath [23] uses entropy power inequalities to solve 

!• certain optimization problems. 

• ^ The entropy power inequality lower bounds the differential entropy of the convolution of 

<S two independent random variables X, Y taking values in M". Its covariance-preserving version 

d states that 

H{VXX + Vl - XY) > XH{X) + (1 - X)H{Y) for all < A < 1 . (1) 

Eq. ([l]) can be shown [221 132] to be equivalent to the more commonly used statement 

^2HiX+Y)/n y ^2H(X)/n _j_ ^2H(Y)/n _ 

The latter explains the terminology as e^^^^)/'^ is the power, i.e., variance of a Gaussian ran- 
dom variable with identical entropy as H(X). Inequalities such as ([I]) are closely related to 
Log-Sobolev inequalities (see e.g., |31|) as well as Brunn-Minkowski-type inequalities |6J. Gen- 
eralizations to free probability [30j and quantum states have been considered. 




Figure 1: The setting of the (unconditional) quantum entropy power inequahty: two indepen- 
dent sets of modes (i.e., a product state px ® pv) combine at a beam-sphtter of transmissiv- 
ityAG[0,l]. 



Among known quantum generahzations is the photon-number-inequahty conjectured in [13] 
(and proved for special cases such as Gaussian states in |12j), and a version of ([I]) involving 
von Neumann entropies S{p) = — tr(plogp) instead of differential (Shannon) entropies [20]. 
Formally, this statement is obtained by substituting states px,pY of n bosonic modes for the 
random variables {X,Y), and using a beamsplitter to process the product state px ® Py, see 
Fig. [I] This results in an output denoted pxmxY on one of the output arms (see below for a 
precise definition), and the corresponding quantum generalization states that 

S{X Ha Y) > XS{X) + (1 - X)S{Y) for all < A < 1 . (2) 

As discussed in [21] , Eq. (J2]) provides strong upper limits on classical communication over addi- 
tive thermal noise channels. Indeed, for suitable parameter regimes, the resulting upper bounds 
on the classical capacity are close to the Holevo-Schumacher- Westmoreland lower bound achiev- 
able by coherent states. This limits the degree of potential additivity violations, implying that 
coding strategies using simple product states are close to optimal for thermal noise channels. 

2 Conditional entropy-power inequalities 

Most applications of the classical entropy power inequality make use of a version for condi- 
tional entropies H{X\Z) = H{XZ) — H{Z). It is easy to see that, if {X,Y) are conditionally 
independent given Z, then 

H{\^X + Vl-XY\Z) > XH{X\Z) + (1 - X)H{Y\Z) for all < A < 1 . (3) 

Indeed, (pi) is an immediate consequence of ([I]) and the fact that the conditional entropy H{X\Z) 
Y^^Pz{z)H{X\Z = z) is the average of the entropies H{X\Z = z) of the conditonal distribu- 
tions {Px\Z=z}z- 

It is natural to ask whether a generalization of (pi) is true for states pxYZ for which the 
n-mode systems X and Y are conditionally independent given the quantum system Z, i.e., 
whether 

S{X Ha Y\Z) > XS{X\Z) + (1 - X)S{Y\Z) . (4) 

We will argue that ^ is useful to estimate entanglement-assisted capacities of additive noise 
channels. The inequality may have additional applications in multi-user quantum information 
theory. 

Establishing an inequality of the form (El) appears to be non-trivial because we cannot 
simply condition on the quantum system Z (unless, of course, it is purely classical). However, 



the following simplification is immediate: it suffices to establishes an inequality of the form 

S{X Ha Y\ZiZ2) > XS{X\Zi) + (1 - X)S{Y\Z2) for all product states pxz^ ® Pyz^ ■ (5) 
This is because any conditionally independent state pxYZ has the Markov form [15J 

j 

and the von Neumann entropy satisfies S{^-PjPj) = J2jPj^iPj) '^H{p) on direct sums, where 
H{p) is the Shannon entropy of the distribution {pj}j- 

Here we prove inequality ^ for all pairs of Gaussian states pxzi and Pyz2- We conjecture 
that Eq. (0) holds in general for arbitrary conditionally independent states px^x^z- If this is the 
case, the implications for entanglement-assisted capacities discussed here extend to all additive 
(but not necesssarily Gaussian) channels. 

It may be possible to find a proof of Eq. Q using a similar strategy as we employ in the 
Gaussian case. However, this will require a novel analysis of the evolution of conditional en- 
tropies under a certain Markovian evolution. In the Gaussian case, this is based on perturbation 
theory for symplectic eigenvalues, as we explain below. 

3 Implications for entanglement-assisted communication 

Quantum communication channels are characterized by different capacities depending on the 
additional auxiliary resources available, as well as whether or not the communicated information 
is classical or quantum. Here we focus on entanglement-assisted classical capacities which 
are arguably best understood. Consider a point-to-point scenario where a sender A tries to 
communicate to a receiver C over a channel E : Bi^A) — )■ B[C\ The entanglement-assisted 
classical capacity Ce{E) is defined as the maximal rate (in bits/channel use) at which classical 
bits can be transmitted reliably if the sender and receiver share an unlimited amount of prior 
entanglement. 

In sharp contrast to the unassisted classical [H] or the quantum capacity |,28j, the quan- 
tity Ce{E) is additive |T]. In particular, it has the single-letter expression Ce{E) = sup^I^S, p) 
in terms of the quantum mutual information 

I{£,p) = S{p) + S{£{p)) - S{{E ® /^O(^AA')) =■■ I{^ ■■ C)(£«/,,)(*..0 ' 

where '^ aa' is a purification of the input density operator p. (This statement is a general- 
ization |21 US] of the Holevo-Schumacher- Westmoreland theorem.) The quantity I{S,p) has a 
number of nice properties: it is positive and concave with respect to the input state p. 

For channels involving infinite-dimensional state spaces, such capacity results can be adapted 
by certain limiting procedures - see [17] for a rigorous analysis of the entanglement-assisted 
capacity. It is necessary to constrain the inputs to obtain meaningful results: one is led to 
consider the quantity C£;(£^,N) obtained by maximizing over all input states p with mean 
photon number upper bounded by N, i.e., 

Ce{S,N)= sup I{S,p). (6) 

p:(ata)p<N 
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Figure 2: The additive noise channel £\,aE with a transmissivity A-beamsphtter and environ- 
ment in the state aE is schematically illustrated by the dotted box. Its entanglement-assisted 
capacity CE{S\,aE^ N) is obtained by maximizing the mutual information I [A : C) over all input 
states pA (with purification '^ aa') subject to the mean photon number constraint tr (a^ap) < N. 



We are interested in estimating ^ when 8 = £x,aE is an additive noise channel. Such a 
channel is characterized by the transmissivity A G [0, 1] and a state aE of the environment, see 
Figure 12} Both the input as well as the environment of the channel £^a,o-b consist of n bosonic 
modes; the two interact with a beamsplitter of transmissivity A. The output of the channel is 
a set of n modes (see Section 8.1 for a precise definition of these expressions). We will focus on 
n = 1 (this often being sufficient because of additivity), although the entropy power inequality 
also applies ioi n > 1. 

The entanglement-assisted capacity of Sx,aE can easily be upper bounded by twice the 



maximum output entropy, i.e.,[^ 



CEiS. 



■^,o-EJ 



N) < 2 max 

p:(ata)p<N 



S{C) < 2(?(N^ax) 



(7) 



Here the second inequality uses the fact [31] that Gaussian states maximize entropy under 
a given photon number constraint: g{x) = {x + 1) log(x + 1) — xlogx is the entropy (in 
nats) of a Gaussian state with mean photon number x, N^ax = AN + (1 — X)Ne is the 
maximal mean photon number at the output, and N^; = {a)a)crj^ is the mean photon number 
of the environment. In the limit A — ?■ 1 of perfect transmission, the rhs. of ([T]) becomes twice 
the entropy of the input as achievable by dense coding [3J. The following Corollary to our 
conditional entropy power inequality (Theorem 8.1 below) improves on the upper bound ([T]). 



Corollary 3.1. Let Sx^^e be the additive noise channel with transmissivity A and environment 
in a state aE with mean photon number Ne- If ^^e is Gaussian, then the entanglement- assisted 
capacity satisfies 



Cs(£:A,.„N)<(7(N^ax) + A^(N) 



X)S{E) 



(8) 



where N^ax = AN + (1 — A)N£;. Moreover, if conjecture Q holds for all states, then the 
bound (|8| holds even in the case where aE is not Gaussian. 



^The proof of (IT]) follows immediately from (see Figure^ 

I{A' : C) = S{C) ~ S{C\A') = S{C) + S{C\E'D) < 2S{C) . 

Here E' purifies the environment, and D is the second beam-splitter output. The second identity uses the purity 
of the overall state on A'CE'D and the inequality is subadditivity of the entropy. 



We emphasize that Corollary |3.1| is of interest mainly in cases where the channel is not 
completely characterized as e.g., in quantum cryptography. Under conjecture (|4]), it gives a 
universal upper bound independent of the detailed structure of the environment's state aE- 

Proof. Assume that aE is Gaussian such that Sx^as is a Gaussian operation. Then the opti- 
mization ^ can be restricted to Gaussian states |18j. Consider an arbitrary Gaussian p^, and 
let "^AA' be a Gaussian purification. Then 

I{A' : C) = S{C) - S{C\A') < S{C) - {XS{A\A') + (1 - X)S{E)) 

by the conditional entropy power inequality for Gaussian states. The claim then follows from 
the maximum entropy principle [3l] (i.e., the fact that Gaussian states maximize entropy under 
a constraint on the second moments) because S'(y4|A') = —S{A) for a pure state "^aa'- 

The case of a non-Gaussian state aE follows in a similar manner (under conjecture (ll])). 
Here we cannot restrict the optimization ^ to Gaussian states. D 

Figure [3] shows a comparison of this bound with the known capacity of the thermal noise 
channel. 

Note that if the state aE is Gaussian, then the channel ^a,ct_b is a Gaussian operation. Var- 
ious capacities of such channels have been studied in detail [IH] (see (TU] for a recent review). 
In particular, the entanglement-assisted capacity C£;(£^A,crB, N) is the result of a convex op- 
timization problem: it is obtained by maximizing the mutual information over the (convex) 
set of covariance matrices M associated with Gaussian states p satisfying the photon number 
constrainlj^ In principle, this can be addressed using efficient numerical algorithms. Further- 
more, in special cases, it is possible to give an explicit expression for C£;(£A,crB,N): Holevo 
and Werner [TS] computed the entanglement-assisted capacity of the attenuation/amplification 
channel with classical noise (the most general one- mode channel not involving squeezing). Simi- 
larly, the entanglement-assisted capacity of the broadband lossy channel was discussed in detail 
in [11]. The crucial feature enabling these calculations is gauge-invariance of Sx^o-e- Because 
the entanglement-assisted capacity is additive, this is equivalent to the state aE being thermal 
with respect to the mode operators defining the beam-splitter. In contrast. Corollary |3.1 gives 



a bound which is applicable to all Gaussian additive channels without further restrictions. 

Recall that the entropy power inequality |2T] provides an additive upper bound on the classi- 
cal capacity of thermal noise channels. Because the entanglement-assisted capacity is additive. 



the conditional entropy power inequality plays a somewhat different role in Corollary |3.1[ it 
substitutes an optimization problem by a bound depending on simple universal parameters 
(i.e., the entropy and the mean photon number of the environment). One may hope that the 
conditional entropy power inequality may also be used to address additivity problems in the 
context of entanglement-assisted communication. A natural candidate problem here is the rate 
region of the quantum multiple access channel (MAC) characterized by Hsieh, Devetak and 
Winter ^\ , where a single-letter formula is not known. The use of conditional entropy power 
inequalities is especially suggestive in the case of the additive bosonic MAC (see e.g., [35]). 
where Alice and Bob are connected to a receiver Charlie via two arms of a beamsplitter. As 
shown by Cezkaj et al. [7j (see also [8]), this scenario exhibits uniquely quantum activation 



■^This follows from the fact that the maximum mutual information in m is achieved on a Gaussian state p 
and the concavity of I{£, p) in p, see [T5] . 





(a) Ne = 1/2, N = 10 and A e [0, 1] 



(b) Ne = 1/2, A = 1/4 and N e [0,4] 





(c) Ne = 0.005, N = 10 and A e [0, 1] 



(d) Ne = 0.005, A = 1/4 and N e [0,4] 



Figure 3: Here we compare the upper bound of Corollary |3.1| with the known capacity of the 
thermal noise channel (with environment photon number N^;). The latter is given by [IH] the 
expressionCi; (gA,a^,N) = ^(N)+^(N^ax)-^( P + N^ax-N-l)/2)-(7((D-N^ax + N-l)/2), 
where D = ^/(N + N^ax + 1)^ — 4AN(N + 1). We are interested in the case where aE is close 
to the vacuum state corresponding to a pure loss channel. We plot the capacity Ce, the EPI 
upper bound ([s]) and the maximum entropy upper bound ([T]). While in this case, the exact value 
of the capacity is known, these figures illustrate how the entropy power inequality improves 
over the trivial bound. Its insensitivity to the exact form of the environment's state may be 
useful in certain applications. In contrast, the expression for Ce depends on the fact that aE 
is a thermal Gaussian state, that is, the gauge-invariance of Sx^^j^. 



effects: providing entanglement-assistance to Bob can boost Alice's maximal rate of communi- 
cation. This is in contrast to analogous classical settings, where providing additional resources 
to Bob cannot change her maximal rate (the latter is determined by her signal power). While 
the conditional entropy power inequality can be used to bound the strength of this activation 
effect, a direct application does unfortunately not appear to yield fundamental new insights 
into the additivity problem for the bosonic MAC. 



Outline 

In Section |4| we recall basic definitions. In Section [5| we perturbatively compute the first-order 
corrections to the symplectic eigenvalues of symmetric matrices relevant in our context. In 



Section[6| we apply these perturbative results to obtain the asymptotic scaling of the conditional 
entropy, as well as its infinitesimal rate of increase when part of a Gaussian state undergoes 
diffusion. In Section [7| we connect this to Fisher information by establishing a de Bruijin-type 
identity for conditional entropies. In Section |8j we complete the proof of the entropy power 
inequality for conditional entropies. 

4 Basic definitions 

Consider A^ bosonic modes described by mode operators R = {Qi, Pi, . . . , Qn, Pn) satisfying 
the canonical commutation relations 

/ 1" ®^ 
[Rk,Re]=iJk,e where -^n = [^ g 

A Gaussian quantum state p = Pj^j ^ on this system is completely specified by its first and 

second moments, i.e., the displacement vector d = {di, . . . ,d2N) ^ I^^^ and its (symmetric) 
covariance matrix M, defined by 

dk = tT{pRk) and Mj^k = tr(p{i?j - dj, Rk - 4}) • 

Here {A, B} = AB + BA. We call pm,o a centered Gaussian state and often write pm for it. 

A Gaussian operation maps Gaussian states to Gaussian states and is determined by its 
action on (M, rf). An example is a displacement (or Weyl) operator W{E,)i S, € M?^: this is a 
unitary operation satisfying 

WiOPM4^iO^ = Pm4H ' (9) 

for all Gaussian states pM,d- It will sometimes be convenient to write the conjugation map as 
WAp) := W{^)pW{^y . Statement ^ is equivalent to the Heisenberg action on the mode 
operators 

W{i)^RkW{i) = Rk + ik for all fc = 1, . . . , 2Ar . (10) 

A matrix S satisyfing SJj\fS^ = J^ is called symplectic. A symplectic matrix S uniquely 
defines a Gaussian unitary Us by the action on mode operators 

2N 

UlR.Us = Y, Sj,kRk =: R'j for all j = 1, . . . , 2Ar . (11) 

fc=i 

Because S is symplectic, the transformed operators R! = {R[, . . . , -R27V) = (Qi? -Pi? • • • ? Q'n^ P'n) 
again satisfy canonical commutation relations. It is convenient to define the (transformed) 
creation and annihilation operators 



1 



The associated number operators 

nk = alak = ^((Ql)^ + (P^)^) -h, /c = 1, . . . , AT (12) 



are mutually commuting, and there is a simultaneous orthonormal eigenbasis of the form \n) 

\ni, . . . ,nk) with nk\n) = nk\n). 



Eq. (11) translates into the action 

UspM,dUl = pSMST,Sd , (13) 

on the covariance matrix and the displacement vector. In particular, displacement operators 
and unitaries of the form Us can be used to diagonalize any Gaussian state (with (|9| and ( 13 )). 



More precisely, Williamson's theorem [33] states that a symmetric positive definite matrix M 
can be diagonalized by a symplectic matrix S, i.e., 

SMS^ = diag(Ai, Ai, As, As, ... , A^, A^v) =: ^iv(A) . (14) 

We call A = (Ai, . . . , Xn) the symplectic eigenvalues of M. Observe that this list may include 
multiplicities (i.e.. A, = \j for i ^ j). We will occasionally denote the set of distinct symplectic 
eigenvalues as Sspec(M). The quantity A(M) = '^^'^xx^zSspedM) xj.x 1-^ ~ -^1 "^i^l be referred to 
as the symplectic gap of M. 



With (13), identity (14) implies that a centered Gaussian state pm with covariance matrix 



M can be brought into product form as 

n —/3h- 

UsPmUI = (g) ^^ , (15) 

where the inverse temperatures Pj = P{Xj) are given in terms of the symplectic eigenvalues by 

/3(A) = log ^. (16) 

Note that /3(A) is monotonically decreasing with increasing A, /3(A) < 2 for A > 2, and 



limA-s.oo /3(A) = 0. If Xj = 1, then the factor ^ J^.i. in (15) needs to be replaced by the pure 
'vacuum' state |0)(0|q',p' associated the the mode operators Q'j, Pj. Eq. (15) shows that the 
number states {|?T-)}„eN^ corresponding to the transformed mode operators R' are an eigenbasis 

of UsPmUI- 

The entropy S{p) = — tr(plogp) of an A^-mode Gaussian state p = pM,d only depends on 
the symplectic eigenvalues A = (Ai, . . . , Xn) of the covariance matrix. To express it, it is useful 
to define the mean photon number N(Afc) in the eigenmode k by the function 

N(A) = (A-l)/2 (17) 

of a symplectic eigenvalue A. Then the entropy is given by 

N 

S(p) = ^^(N(Aj)) with ^(N) := (N + l)log(N + l)-NlogN . (18) 

Let us give a simple bound on the dependence of the entropy on the symplectic eigenvalues. 
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Lemma 4.1. Suppose p and a are N-mode Gaussian states with (arbitrarily ordered) symplectic 
eigenvalues u = (z/i, . . . , z/tv) and X = (Ai, . . . , \n), respectively, and assume that \j ^ 1 for all 
j = 1, . . . ,N. Let ||i7 — A||oo = niaxi<j<7v \i^j — Xj\ and A^, = mmi<j<Ar Xj. Then 

\S{p) - S{a)\ < I M|z7- A|U/3(A.) + ^^^ 

The requirement that a has no symplectic eigenvalue equal to 1 (or equivalently A* > 1) implies 
that none of the eigenmodes factors out in a pure product state, and we can usually assume 
this without loss of generality in our considerations below. 



Proof. We have g'(N) = log((N + 1)/N), and combining this with (17), we obtain the inverse 
temperature according to 

^'(N(A)) = /3(A) . (19) 

Furthermore, we have 

/3'(A) = -3^. (20) 



Let us define the function G(A) := c/(N(A)). Because N'(A) = A/2, Eqs. ^ and ([20]) imply 

G"(A) = /3(A)/2 and ^'(A) = -^^ • 

Hence the Taylor series expansion gives 



|G(A + e)-G(A)|<e/3(A)/2 



e2 



2|A2-1| ■ 
Since (3 is monotonically decreasing, we get 

1 ^ \ n/ ■ \ \ niaXfc I Z/fc — A 



max|G(z/j) - ^(Aj)! < -max|z/fc - Afc|/3(minAj) + . „ n 
j Z k j zmm,' A,- — i 

The claim follows from this because S{p) = J2j=i ^i^j) ^^^ S{a) = J2j=i ^i^j) ■ ^ 

5 Perturbation theory for symplectic eigenvalues 

We begin with a straightforward application of degenerate perturbation theory which is illus- 
trative of the method. We will need slightly more involved statements for bipartite systems 
below (Lemma [5.2 ). Note that a similar perturbative analysis was used in a different context 
in [26l Appendix B]. 

Lemma 5.1 (Perturbation of symplectic spectrum to 0th order). Let M be a covariance matrix 
of N modes and consider the symplectic eigenvalues A(e) = (Ai(e), . . . , Xn{^)) of 

M°°(e) = I + eM , 

where we assume that \\M\\ = 0(1) and e ^ 1. Then 

Aj(e) = l + 0(e) forallj = l,...,N. 



Proof. The symplectic eigenvalues A(e) can be obtained from the spectrum spec(zJjvM°°(e)) of 
iJiqM°^{e) since this operator has eigenvalues (Ai(e), — Ai(e), . . . , Aiv(e), — AAr(e)). Let us write 

iJNM'^ie) = H + V where H = Un and V = ieJ^M . 

Observe that the eigenvalues of H are {1,-1} and ||V^|| = 0(e). This means that we can apply 
standard degenerate perturbation theory: the spectrum spec(if + V) is given by 

{/i + spec{V\nJ I /i G spec{H)} + 0{e') (21) 

to first order in e, where ^|^^ denotes the restriction of V to the degenerate eigenspace 7{^ of H 
to eigenvalue fi. An immediate consequence is that (by assumption on ||M||), the symplectic 
eigenvalues of M°°{e) are of the form Aj(e) = 1 + 0(e) for all j = 1, . . . , A^. D 



To obtain the correct constant (i.e., the first-order correction in (21)), it is necessary to 



compute the restriction V\-Hf, to the eigenspace with eigenvalue /i. For this purpose, it is 



convenient to use a basis of Ti^. For example, for the case of Lemma 5.1, we can define the 
vectors 

K) = ^(^, 1) = ^m + |2)) and \v-) = i=(l, t) = i=(|l) + t\2)) 

in M^ (Here and below, we will often use \j) to denote standard orthonormal basis vectors 
in M^-'^.) Importantly, these vectors satisfy 

iJi\v^) = ±\v^) , (22) 

and 

{v^\iJiDi{a)\v^) = ±a {v'^\iJiDi{a)\v^) = for all a G M . (23) 

Writing (0, 0)®'' = (0, 0) © ■ ■ ■ © (0, 0) {k summands), we can define 



\v 



± 

'j,N 



) = (0,0)(^'-^)©|t;±)©(0,0)®(^-^'-^) . 



where \v ) is at mode j (of A^ modes). 



TV- 



Because of (22), the eigenspace 7ii of if to eigenvalue 1 has orthonormal basis {\v^j^)}j=i 
Similarly, the eigenspace "H-i has orthonormal basis {\v~n)}j=i,...,n- The restriction of V to "Hi 
is given by the matrix {{v~l'\V\v^))j^k=i,...,N, and diagonalizing this matrix gives the first-order 
corrections to the eigenvalue 1 (the degeneracy is lifted). We will omit a more detailed discus- 
sion here as we will need a more general version (including a reference system B). However, the 
proof of the following Lemma proceeds in this fashion; the key here is to apply perturbation 
theory to a suitably transformed matrix. 

Lemma 5.2 (Perturbation of symplectic spectrum to 1st order). Let 

he the covariance matrix of a state on m + n modes with respect to the ordering 

^ = [Ql 5 -'1 } ■ ■ ■ } Qmi ^mi Ql 1 ^l ■,■■■■, Qn ; -< n ) 

of modes. 
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(i) Assume that \\Mab\\ = 0(1) and e ^ 1. Then the symplectic eigenvalues A(e) 
(Ai(e), . . . , Am+n(e)) of the covariance matrix 

MTsi^) = ( ^^ °') + eMAB =■■ Ia®Ob + ^Mab (24) 



^0 

are of the form 

\{^) = Vj{^) + O(e^) for j = l,...,m and Xm+j{e) = euj + O(e^) for j = l,...,n 

where f]{e) = (rji^e), . . . , '^.m(e)) are the eigenvalues of I a + ^Ma and U = (i/i, . . . , z/„) are 
the eigenvalues of Mb- 

(a) Let X = (Ai, . . . , Xm+n) be the symplectic eigenvalues of Mab o^nd let Sab be the symplectic 
matrix diagonalizing Mab, ^-e., 

SabMabSab = Dm+n{X) . 

Suppose e <^ A{Mab)- Let A(e) = (Ai(e), . . . , Am+„(e)) be the symplectic eigenvalues of 

M%{e) = MAB + elA®QB ■ (25) 

Then degeneracies are split according to 

\{i\X^ = X}\ = \{e\ Xe{e) G [A - A/2, A + A/2]} (26) 

where i ranges over i E {1, . . . ,m + n}. Furthermore 

J2 He)= J2 {^+l^4SAB{lA®0B)S'^Bf^)+O{e') 

£:A£(e)e[A-A/2,A+A/2] £:A£=A 

for each X G Sspec{M ab) ■ Here [Z]^^'^ denotes the submatrix corresponding to the i-th 
mode, i.e., the 2 x 2-matrix with entries {[Z]^^^)ij = Z2i-i+i,2e-i+j , i,j ^ {0, 1}. 



Proof of Lemma 5.2 (|ii: 

Our goal is to compute the symplectic eigenvalues Ai(e), . . . , Am+n(e) of the covariance matrix 

M^si^)=lA(BOB + eMAB , 

to first order in e. Let i? = (i/i, . . . , z/^) be the symplectic eigenvalues of Mb, and let Sb be the 
symplectic matrix diagonalizing Mb, i.e., 

SbMbSI = D„(z7) , 

where Dn{v) = diag(i^i, z^i) © ■ ■ ■ © diag(z/„, z/„). Similarly, let a(e) = (ai(e), . . . , arn(e)) be such 
that the symplectic eigenvalues of Ia + ^Ma are equal to (1 + ai(e), . . . , 1 + ^^(e)). Observe 
that 

|a,(e)|=0(e) for j = 1, . . . ,m . (27) 
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Let Sa{^) be the symplectic matrix diagonalizing Ia + eM^, i.e., 

SA{e){lA + eMA)SA{ef = Ia + D^{a{e)) . 
Finally, let Sab{.^) = Sa{^) ® Sb- The symplectic spectrum of M^(e) is identical to that of 



MTsi^) := SABie)MUe)SABier = { /^rT c f.^T 



T _ (Ia + Dm{a{e)) tSA{e)LABSB 
eSsL^BSAitf eDrr{iy) 



In particular, the eigenvalues of the matrix iJm+nM'^^{e) are given by 

spec{iJm+nM^Bi^)) = (-^i(e), -Ai(e), . . . , A„+„(e), -A„+„(e)) . (2^ 

Since Jm+n = Jm ® Jn, '^^ obtaiu 

aMr.(e) = H^V where H = ('^™ °) and V = (■'^•}5f*>' ^^]^. 

where we introduced the abbreviation X(e) = ieSA{€)LABS'^- By the assumption ||M^b|| = 
0(1), we have ||X(e)|| = 0(e) and \\Dn{v)\\ = 0(1). Similarly, we have ||D^(a(e))|| = 0(e) 



because of Eq. (27). We conclude that ||V"|| = 0(e). On the other hand, the eigenvalues of H 
are {0, 1, —1}, hence we can apply first-order (degenerate) perturbation theory to compute the 
spectrum of iJM'^^{e). 

We consider the eigenspaces separately. 

• /i = 0: Consider the eigenspace "Ho of H to eigenvalue fj, = 0. It is easy to check that the 
vectors {t;^^, v^j^, . . . , v^^, v'^b) defined by 

HB)-=OA®\vfj 

are an orthonormal basis of "Hq- (Here we write 0^ for the zero- vector (0, 0)®™.) Since 
we get the matrix elements 

{'Vj,B\yK,B) = {vl^\itJnDn{v)\vl^n) = ^j^k^aA^ ' ^^j) 



according to (23). In particular, the restriction \^|-Ho is described by a diagonal matrix 



with diagonal elements of the form {±ez/j}"^^. With (21), we have found m symplectic 
eigenvalues of the form 

Xm+j{e) = euj + 0{e^) for j = 1, . . . , n . (29) 

(Only the non-negative entries on the diagonal are symplectic eigenvalues according 
to ([28]).) Note that these are simply the eigenvalues of Mb, up to a factor of e. 
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/i = 1: An orthonormal basis of Hi is given by the vectors {|f^^)}"l^ defined by 



4a)-=HJ®^ 



B 



where 0^ stands for (0, 0)®"'. It is easy to check that 

i.e., the restriction of V to "Hi is diagonal when expressed in this basis. In conclusion, we 
found m symplectic eigenvalues of the form 

Aj(e) = 1 + aj{e) + 0{e'^) for j = 1, . . . , m . (30) 

By definition of aj(e), these are equal to the symplectic eigenvalues rij{e) = 1 + aj(e) of 
I A + ^Ma in order 0(e). 

u = — 1: Here we find an orthonormal basis with vectors IvfA) := \vf) © Or. The 
restriction of V to "H-i is diagonal with diagonal entries (— ai(e), . . . , — am(e))- The 



corresponding eigenvalues are the negatives of (30), consistent with (28) 



Proof of Lemma 5.2 ([n]) 



Let us briefly recall the definitions involved in the statement. We consider the covariance matrix 



I A 





M%{e) = Mab + e n n = ^ab + tU © 0^ , 



where Mab has symplectic eigenvalues A = (Ai, . . . , Xm+n) and is diagonalized by Sab, i-e., 

SabMabSab = Dm+n{X) ■ 

We assume that e ^ A, where A = /S.{Mab) is the symplectic gap. Then M^^(e) has the same 
symplectic eigenvalues as 

M%{e) = D^+niX) + tSAB{lA ® Qb)SIs . 

In particular, it suffices to find the eigenvalues of 

iJm+nM^^{e) = iJm+nDm+n{X) + iedm+nSABi^A © Ob)Sj^q =\ H + V . 

Because H := iJm+nDm+nW has gap A, ||^Jm+n'S'^B(-^A©0_B)5'J^|| = 0(1) and our assumption 
e ^ A, we can apply (degenerate) perturbation theory to compute the spectrum of zJM^^(e). 
The operator H has eigenvectors {wf, . . . , w^+„) = {v^a^ • • • ; "^m A; '^i'b; • • • ? "^n b) "with 
eigenvalues Ai, . . . , Xm+n- In particular, for A G Sspec(MAB), the eigenspace "Ha of H is 

Tix = span{|u7^) | 1 < i < m + n with Xi = X} . 

Suppose that restriction V\-}i^ of V := ieJm+nSAB{lA®OB)SAB ^o ^a has eigenvalues 6i, . . . , OdimUx- 
According to (degenerate) first-order perturbation theory, the matrix iJm+nMj^^^e) has dim?^A 
eigenvalues of the form 

A + % + O(e') for J = 1, ... , dimT/A , (31) 
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where 6j = 0(e). Furthermore, the hst (31) includes all eigenvalues Xi{e) of iJm+nM^g{e) in 



the interval [A — A/2, A + A/2] and the number of such eigenvalues is equal to (26) by the 
assumption e ^ A. We conclude that 



dim "Hi 



J2 A,(e) = (dim?^;,) ■ A + ^ 6,+ 0{e') 

^:Af(e)e[A-A/2,A+A/2] j=l 

= (dimnx) ■ A + tr(K|«J + 0(e2) (32) 

Because the restriction V\-u can be expressed as 

j,k:\j=Xk=X 

we obtain 

tr{V\nJ =ieY, {wj\Jm+nSAB{lA © Qb)SIbH) . (33) 

£:A,=A 



Combining (p33|) with dim'HA = |{^ ^ A^ = A}| and Eq. (32) gives 



Y, \,{e)= Y, {\ + ie{wi\J^+nSAB{lA®QB)SlB\wt))+0{e') . (34) 

£:A^{e)e[A-A/2,A+A/2] t\i=\ 

Since \wj) = ^(«|2£ - 1) + \2i)) and J^+„ = ^;i+" \2k - l){2k\ - \2k){2k - 1|, we have 
{w^\Jm+n = 75 (^(2^*1 + {2i — 1|), and this takes the form 

{w^ \Jm+nSAB{lA © Ob)5'^b|'U;^ ) = - (^^2^,2£-l — "^^2^-1,2^ + ■^2£,2^ + ■Z'2^-l,2^-l) 

where Z = Sab{Ia © ^b)S^b. Since Z^ = Z is symmetric, this is equal to 

{wi\Jm+nSAB{lA © 0^)51^10 = ^ ^[^^^(/a © Ob)^!^]^') , (35) 



and the claim follows by inserting ( 35 ) into ( 34 ) . 



6 Conditional entropy and diffusion 

A central tool in the proof of the quantum entropy power inequality [20] is the one-parameter 
semigroup {e^^}t>o of completely positive trace-preserving maps generated by the 'diffusion' 
Liouvillean £ (see [20] for a definition of the latter). For all t > 0, the map e*^ is Gaussian and 
can therefore be defined in terms of its action on the covariance matrix M and the displace- 
ment vector d: A Gaussian state p described by (M, d) is transformed into a Gaussian state 
p{t) = e*^(p) with covariance matrix (M(t) = M + tl,d), i.e., the transformation governing the 
evolution for time t is 

p !-)• e*^(p) 

^ ^tc ^ for all covariance matrices M and displacement vectors d . 

{M,d) U {M + tI,d) 
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In this section, we revisit and extend statements of [20j about the behavior of the entropy S'(e*^(p)) 
as a function of time t. In particular, we specialize to Gaussian initial states p and extend our 
considerations to conditional entropies. 

More precisely, we consider the case where diffusion acts only on a subset of modes. Con- 
cretely, assume that our system is bipartite, with system A consisting of m modes, and system B 
consisting of n modes. Diffusion for time t acting on the modes in A only is described by the 
superoperator e*"^^ ® Ib (where Ib is the identity superoperator on B). In particular, this 
family of superoperators is specified by the transformation 



PAB ^ (e*^^ ®/b)(pab) =:pAB(t) 

(MabJab) "'''^'^ {MAB + t{lA(BOB)jAB) 



(36) 



for all covariance matrices Mab G Mat2(m+n)(I^) and displacement vectors (Iab G M2('"+"). We 
will examine the conditional entropy S{A\B)p^gi^t) for the evolved state pABit), given some 
Gaussian initial state pab = Pa_b(0). In Section 6.1, we show that S{A\B)p^g(^t) scales as 



we 



a universal function for t — t- oo (independent of the initial state Pab)- In Section 6.2 
derive an explicit expression for the infinitesimal rate of change of this quantity in terms of the 
covariance matrix of pab- 

6.1 Scaling of the conditional entropy in the infinite- time limit 

In the infinite-time-limit, the entropy of the time-evolved state e*^(p) scales as a universal 
function of time t which is independent of the initial state p. This statement was shown for 
general states in jJOl Corollary 3.4]; here we give a simple argument for Gaussian states for 



completeness. We will also need this statement in the proof of Lemma |6.2| which deals with 
conditional entropies. 

Lemma 6.1 (Scaling of entropy in the infinite-time limit under diffusion). Let p be an (arbi- 
trary) Gaussian state of N modes. Then 

\im\S{e'^{p))-N.g{{t-l)/2)\=0. 

Proof. Let M be the covariance matrix of p. The covariance matrix of p{t) = e*^(p) has the form 
M + tl = tM°°{l/t). Therefore, the symplectic eigenvalues are of the form Vj = t{l + 0{l/t)) = 



t + 0(1) according to Lemma 5.1 For t > 1, the matrix tl is a valid covariance matrix with 
symplectic eigenvalues Xj = t ioi j = 1, . . . ,N. Let a{t) be the centered Gaussian state with 



covariance matrix tl. Lemma 4.1 applied to p(t) and a{t) gives 



\S{p{t)) - S{am < |0(1) (/3(t) + ^^ ^0 for t ^ oo . 

Since S{a{t)) = N ■ ^(N(t)) = N ■ g{{t - l)/2), the claim follows. D 

A similar statement holds for conditional entropies: 

Lemma 6.2 (Scaling of conditional entropy in the infinite-time limit under diffusion). Let pab 
be a Gaussian state of {m + n) modes A and B . Define pABif) = (e*^-* ® Ib){.Pab)- Then 



t—^oo 



\im\S{A\B)p^^^,)-m-giit-l)/2)\=0. 
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Comparing Lemma 6.2 with Lemma 6.1 suggests that (e*^-^ ® Ib){pab) approaches a product 
state for large times. Quantifying this convergence for general (possibly non-Gaussian) states 
may provide an avenue to proving conjecture |4} 

Proof. Let u = (z/i, . . . , i/) be the symplectic eigenvalues of the covariance matrix Mb- Note 
that we can assume without loss of generality that 

Uj 7^ 1 for all j = 1, . . . , n . (37) 

Indeed, if we had Uj = 1 for some j G {1, . . . , n}, then the corresponding eigenmode in B is in 
a pure state and the state pab factorizes. This is true also for the time-evolved state pABif) 
(because B is unaffected by the evolution), hence such eigenmodes do not contribute to the 
entropy S{A\B)p(t) and can be traced out. 

The covariance matrix of the state pAB{t) can be written in the form 

Mab + tiA ®0b = til a ®0b + l/t ■ Mab) =: tM^Bi^/t) . 



Applying Lemma 5.2 m with e = 1/t (and multiplying the resulting eigenvalues by t), we 



conclude that its symplectic eigenvalues are 

Xj = Tjjit) + 0{\/t) for j = 1, . . . , m and \m+j = Vj + 0(l/t) for j = 1, . . . , n , 

where {?7j(t)} are the symplectic eigenvalues of Ha + Ma- This means that up to order 0(l/t), 
the symplectic spectrum of pAsit) is identical to the one associated with the product state 
e*^(pyi) ® Pb- In summary, we have the list of symplectic eigenvalues 

pABJt) e*^^(PA) I PB 

m) + o(i/t), u{t) + o(i/t)) ffit) v{t) 

V V 

=:fi =:i> 

Let Pa be a Gaussian state of m modes with symplectic spectrum f], and let pB be a Gaussian 
state of n modes with symplectic spectrum z/. Then 

S{pAB{t)) = S{pa) + S{pb) - (38) 



We can apply Lemma 4.1 to the states pa and e*^'*(p^). Since the latter has symplectic 



eigenvalues r/(t) = (?7i(t), . . . ,?7„(t)), we get 

\S{pa) - S{e"'{pA))\ < y (0(l/t)/3(r/,(t)) + 0(l/t^) ^^^^^^, _ J ^0 as t ^ oo . 

Here we used the fact that t]j{t) = t + 0(1) for all 1 < j < m, hence rj*{t) = miuj rjjit) — )■ oo for 



t — )■ oo. Similarly, applying Lemma 4.1 to the states pb and pB (and using assumption (37)) 
gives 

\SiPB) - S{pb)\ < ^ (oil/t) ■ /3(z/,) + 0(1/^) . ,\ . ] ^0 for t ^ oo . 

2 \ mmj {Uj — 1 1 / 



Inserting these upper bounds into (38) and using the triangle inequality then gives 



|5((e*^-^ ® Ib){pab)) - Sie'^'^ipA) Pb)\ ^ for t ^ oo 
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Because the states pAsit) = {.^^'~'^®Ib){.Pab) and aABif) = ^^'^^{.Pa)®Pb have the same reduced 
density operator pBit) = o^sit) = Pb, the previous statement imphes that the difference between 
their conditional entropies also vanishes in the limit, that is, 

\S{A\B),^^^t) - S{A\B)^^t)\ ^0 for t ^ cx) . (39) 

Since (TAB{t) is a product state, we have 

SiA\B)^^t) = SiaAit)) . 

The claim then follows from the triangle inequality, i.e., 

|^(A|S),,,(,) - m . g{{t - 1)/2)| < |5(A|5),,,(,) - SK(t))| + \S{aA{t)) - m ■ g{{t - 1)/2)| 



(40) 



because the first term on the rhs. goes to for t — > oo according to (39) and (40), whereas the 
second term goes to according to Lemma 6.1 , D 



6.2 Rate of increase of the conditional entropy 

Next we compute the infinitesimal rate of increase of the conditional entropy under the pro- 



cess (36). 



Lemma 6.3 (Rate of conditional entropy increase under diffusion). Consider bipartite sys- 
tem AB of m + n modes. Let pab be a Gaussian quantum state whose covariance matrix 
Mab has symplectic eigenvalues A = (Ai, . . . , Xm+n)- Let Sab be a symplectic matrix such that 
SabMabSIb = Dm+nW- Define pABif) = {e''^^ ® Ib){pab)- Then 



d 
dt 



m+n 



^J{A\B),,,^,) ^ 



\Y.P{\dtT{[SAB{lA®QB)SlBf^) 



where \Z^^^ is the 2x2 submatrix corresponding to the i-th mode, ([ZY)ij = Z2e~i+i,2e-i+j for 

z,ie{0,i}. 

Proof. Because the reduced density operator tr^ Pab (i) = Pb is independent of time, we have 

d 



d 
dt 



..o^'-^l^''-'" dt 



t=o 



SiAB),,,^,) . 



(41) 



To evaluate the rate of change of S{AB)p^g(^t), let (Ai(t), . . . , \m+n{t)) be the symplectic eigen- 
values of pAB{t). According to expression (|l8| for the entropy, we have 



m+n 



s{pAB{t)) = Y,9m>^m)= E E amum m 

i=l A6Sspec(A/Ai3)^:A^(t)e[A-f ,A+f] 

where we reexpressed the summation using the symplectic gap A = /\{Mab)- With the ex- 



pression (17) for the mean photon number, Eq. (19) and the chain rule for differentiation, we 
have 



9(N(A(«))=9(N(A(0))) + 



^^^.A'(0) + O(*^) 
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Observe that pAsit) has covariance matrix M^^{t) and we can restrict our attention to times t <^ 
A(Mab) without loss of generahty. Hence we can apply Lemma [5.2| ([n]). We obtain 

J2 gmMm = ^(N(A)) ■ \{e \Xe = A}| + ^ 5^ tr[^^^(/^ © 0^)^^]^^) + 0{t') 

£:AHt)e[A-|,A+f] «:A,=A 



for any A G Sspec{MAB)- Taking the derivative at t = and inserting into (42) therefore gives 
d 



dt 



S {pAB{t)) = \ Y. /^(^) E ^^\^ab{Ia © ^b)S\j,, 



T iW 

t=0 4 

AGSspec(A/As) ^;A^=A 



which is the claim because of (41). D 



7 Diffusion, translations and Fisher information 

A key element in the proof of the classical entropy power inequality is de Bruijin's identity; 
it relates the infinitesimal rate of entropy increase to the Fisher information of a family of 
translated distributions. In [20], a quantum version of this statement in terms of the diffu- 
sion semigroup and phase space translations was given. Here we derive a generalization of 
this statement for conditional entropies (but specialized to Gaussian states). Our proof pro- 
ceeds by direct calculation and does not involve any technical subtleties associated with formal 
computations involving infinite-dimensional systems. 

We begin by recalling the relevant definitions. Consider a one-parameter family {p^^-'jeeiR 
of states depending smoothly on the parameter Q. The divergence-based Fisher information of 
this family (at Qq G M) is defined as the quantity 

J(p(^);e)|,=,„ = ^^(pWl|p(^))|,=,„, 

where S'(p||cr) = tr(plogp — plogcr) is the relative entropy or divergence. Two straightforward 
but important consequences of this definition are the reparametrization identities 

J{p^''^^ d)\e=^ = c'j{p^'^- e)\e=o and J(p(^+^); d)\e=e, = J(p('+^^ e)\e=e,+c (43) 

for c G M and its additivity 

J{pf ® pf: 0)\e=e, = ApT; e)\e=eo + J{pf: e)\e=e, ■ (44) 

The latter follows from the additivity of the relative entropy under tensor products. Further- 
more, because of the monotonicity S{£{p)\\£{a)) < S{p\\a) of the relative entropy under CPTP 
maps £ the Fisher information also satisfies monotonicity (see [20]), i-e., 

J(^(p(^));e)|e=o<J(p('^^)k=o. (45) 

de Bruijin's identity involves the family {p^^'-^^fe) = We|A:)(p)}eeK of states obtained by trans- 
lating an A^-mode state p in the direction Rk in phase space, where A; G {!,... ,2iV}. Let us 
write 

27V 

J{p):=Y.J{p'^''''^^Me=o (46) 



fc=i 
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for the sum of the corresponding Fisher informations, de Bruijin's identity (shown in [20] ) 
relates this to the rate of entropy increase under diffusion, i.e., 



d 

dt 



1 



t=o 4 



(47) 



In this section, we derive a version of (47) for Gaussian states which involves an auxiliary system: 



it quantifies the rate of increase in the conditional entropy S{A\B) when A undergoes diffusion. 
In contrast to [20], the proof given here proceeds by direct computation. In Section 7.1, we 



compute the Fisher information of a family of states obtained by translating a Gaussian in 



phase space. In Section 7.2, we combine this with Lemma 6.3 to prove the de Bruijin identity. 



7.1 Conditional Fisher information of translated Gaussian states 

Let pj^j- ^ denote an A^-mode Gaussian state with covariance matrix M and displacement d G 

M^^. Suppose S" is a symplectic matrix such that SMS'^ = D{'y) is diagonal. Let 6 G M^^ be 
arbitrary. We will need the following formula for the relative entropy of p^^ j- and a displaced 

state Pjv/d+e- ^^ have 



1 ^ 



(4J 



i=i 



where F is a function of the symplectic eigenvalues only. 



Proof. By the invariance 5'(f/pf/"''||f/crf/"'') = S{p\\a) of the relative entropy under unitaries, and 
applying displacement operators as well as the unitary Us, we have 

SipM,d\\pAI,d+e) = SipMfi\\pM,e) = SipM-9\\pM,o) = Sipj^(^^)_sg\\pD{y)fi) ■ 

It hence suffices to analyze S{p£){'y),ff\\pD(-y),o), where ff = —S9 G M^^. Because 

PDif),ff = PDi(7i),(r,i,r,2) ® " " " ® PDi{'yM),iV2N-l,V2N) 

is a product state, we obtain (using the additivity of the relative entropy for product states) 



N 



S{PM,d\\PM. 



d+0) 



S{PD{i),fj\\PD{^),o) -2_^S{pD^{^^),{r,2j.un2j)\\PDil,^3)fi) ■ 



The claim therefore follows from Lemma 17.11 



D 



Lemma 7.1. Let ff = {f]Q,f]p) G M^,7 G M, and let D = diag{'y,'~f) be the covariance matrix of 
a single mode Gaussian state. Then 



SipD,rf\\pD,o) = 2 (N + vl + Vl) - ^(N) - log(l - e- 



(49) 



where the mean photon number N and the inverse temperature j3 are given by (17) and (16), 
respectively. 
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Proof. For brevity, let us write pu = Pd,o foi' the centered state. Then we have 

SipD,n\\pD) = -S{pd) - tl{pD,ff\ogpD,o) = -5'(N(7)) - tT{pD,ff^OgpD,o) ■ (50) 

To compute the latter term, we use the expression pc = ^^i^-pA) = (1 ~ e~^)e~^'^, where 
h = a^a = ^(Q^ + P^ — 1) is the number operator and (3 = (3{'y) = log(7 + l)/(7 — 1) the inverse 
temperature. By taking the logarithm, one gets 

- tr (pd,^ log pd) = -log(l -e"^) + /3tr(pD,^) . 

Using the fact that pD,ff = W{fj)pDW{f]y for the Weyl operator W{f]) and the fact that 



according to (10) and Definition (12), we get 



- tT{pD,rf'^ogpD) = - log(l - e-^) + ^ tr(pz)((Q + ^Qf + (P + ^pf - 1/2)) 

/3 



log(l - e-^) + - (tr(pBn) + r/ J + rfp) . 



In the last line, we made use of the fact that po is centered. The claim follows by combining 
this with ([50|. D 



With (48), we can easily compute the Fisher information of a family of displaced states. 



Lemma 7.2 (Fisher information of displaced states). Let Mab be the covariance matrix of a 
centered state pm,o ofm + n modes, where SabMabS^b = Dm+n{l)- Fix some A; G {1, . . . ,2m} 
and consider the family of states {p'-^'^'^-'jeeK, 



obtained by displacing the state Pm,o ^^ the direction Rj. by an amount 6 eM.. Then 

ra+n 



(51) 



J{p 



{0,Rk). 






Proof With (48), we obtain 



„2 m+n 

'5(pM,o||p(''^'^)) = F(7) + -Y1 /3(7,)(^2Vi,^ + ^: 



2j,k) ■ 



The Fisher information is the second derivative of this quantity with respect to at ^ = 0, 
hence the claim follows. D 
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7.2 The de Bruijin identity for conditional entropies of Gaussian 
states 



It will be convenient to define the conditional Fisher information 

2m 

fc=i 



(52) 



by summing over the modes corresponding to system A only. Observe that many properties of 
the Fisher information carry over to this definition: for example, we have monotonicity 



J{£{A)\B) < J{A\B) 



(53) 



for any CPTPM acting on £^ where these quantities are evaluated on the states pab and 
{8 (g) Ib){pab), respectively. 



By combining Lemma 6.3 and Lemma 7.2, we obtain a proof of the following statement. 



Theorem 7.3 (de Bruijin identity for Gaussian states and conditional entropy). Let p^s be a 
centered Gaussian state of m + n modes. Define pABif) = (e*^'* ® Ib){pab)- Then 



d 



^J{A\B),^,,^,) = -J{A\B),,, . 

Proof. It is straightforward to check that S{Ia © 0^)5*"^ has diagonal elements of the form 
(SilA © OB)S'')e,e = ES Slk- Hence 



d 
Jt 



m+n 



S ((e*^-^ ® Ib){pab)) = \Y. ^(^i) (("^(^^ ® Ob)5^)2,-i,2,-i + {S{Ia © Os)5'^)2,,2,) 

^ m+n 2m 



j=i k=i 



according to Lemma 6.3 We conclude from Lemma 7.2 and Definition (52) that 

d_ 
dt t=o 



^S {{e'^- ^ IbKpab)) = ^J{A\B) 



PAB 



The claim then follows because pB{t) = Pb does not depend on time (cf. ([4l|). 
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8 The entropy power inequality for conditional entropy 

Having established the de Bruijin identity for conditional entropies as well as the asymptotic 
scaling of the conditional entropies under diffusion, it is straightforward to prove the entropy 
power inequality for conditional entropies. Indeed, this follows the pattern of known classical 
proofs [9], with minor modifications because we are considering conditional entropies. It relies 
heavily on the Fisher information inequality (a consequence of data processing, as shown by 
Zamir [36]). 



We introduce the necessary definitions in Section 8.1 The conditional Fisher information 



inequality and the conditional entropy power inequality are derived subsequently in Sections |8. 2 
audio 
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Figure 4: This quantum circuit circuit defines the states px^mxX2 ^'^'^ P(x^mxX2)E^E2 fo^ ^^Y 
Droduct input oy, k, ^ Ox^r^. 



product input pxiEi ® pXiEi- 

8.1 Beam splitters, product states and auxiliary systems 

Consider two systems Xj, j = 1,2 with A^ modes each and associated mode operators {Qj^ , Pk }k=i- 
A beam-sphtter with transmissivity A G [0, 1] acting on X1X2 is the Gaussian unitary Us^ de- 
scribed by the symplectic matrix 

S^= n ^T r\T (54) 



Vl^XhN -VXL 



2N 



with respect to the ordering {Q[^\ Pf \ ■■■,Qn\ Pn\ Qf\P?\ ■■■,Qn\ Pjv^) of modes. We 
are interested in the beam-sphtter map £\ = S^^^^^"^ , which is obtained by letting Xi,X2 
interact according to Us^, and discarding the second set of N modes. That is, it is a map from 
2N input modes to A^ output modes; we call the latter Y. Formally, the map Sx is defined as 



^a(pxiXi) = trx2 [Uxpx^X2Ui 



where tr^a denotes the partial trace of over the second set X2 of modes. We will denote 
the output system (i.e., the set of modes Xi at the end of this process) by Y, i.e., think of 
Sx = S^^ '^'^ as a map from systems X1X2 to an output system Y (of A^ modes). Since the 
partial trace tr^a is a Gaussian map, the map £x is Gaussian and completely determined by its 
action on covariance matrices and displacement vectors. This is 

PX1X2 '->^a(PXiX2) 



where 



My = \Mx, + {l-\)Mx2 + ^W^{Lx,x2 + Ll^x2) 
dy = yXdxi + Vl — \dx2 ■ 



(55) 



This follows immediately from (54). 



We will consider input states that are products (across the bipartition Xi : X2). It will be 
convenient to introduce the following maps and states as summarized in Figure |4j 
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the state pximxX2^ given px^, Px2'- Fo^' product inputs pxi ® PX21 the transformation (55) 
specializes to 

Px, ^PX2 ^ (Ty:= S^'^'^^'^ipx, ® PxJ 

where we assume that pxj is a Gaussian state described by {Mx^, dxj), j = 1, 2. We wiU 
denote the output state ay obtained in this fashion by pximxX2- 

the state P{XimxX2)E-i_E2 given pxiEi and px2E2'- More generaUy, consider the map IeiE2 ® 
£^^ ^^ , where Ej^ j = 1,2 are auxiliary systems with L modes each. We will only need 
a description of its action on product states pxiEi ® Px2E2^ where px e are Gaussian 
states with covariance matrix and displacement vectors 

Mx,E, = i^'"' ^^""'^ d, = [dx.jE,) for J = 1, 2 . (56) 

V ^]Ej Ej J 

It is straightforward to verify that this is given by 

PXiEi ® PX2E2 ^ (^YEiE2 ■= {IeiE2 ® ^x'^''^^){pXiEi ® PX2E2) 
{Mx^Er © MX2E2, {di, d2)) ^-^ {MyEiE2, {dy, rf^i, rfijj) 

where 



/XMx, + {l-X)Mx2 V\Lx,E, VI - \Lx2E2 
Mye,E2 = ^Ll^E, Me, 

V Vl^\Ll^E2^ ^ Me2 

dy — y/Xdxi + -\/l — Xdx2 ■ 



(57) 



8.2 The conditional Fisher information inequaUty for beamspUtters 

In [20l Lemmas 3.2 and 6.1], it was shown that the beam-splitter map is compatible with both 
diffusion and translations in the following sense. For all t > 0, wi, U72 G M, G M^™', we have 
the following identitities of Gaussian maps: 

gX,X2^Y ^ ^^tc^^ ^ ^tc^^ ) ^ ^tcy ^ ^x,X2 (58) 

gX,X2^Y ^ /yyXi^ y^X2^ ^ y^Y Q gX,X2^Y ^j^g^g ^ _ ^ Vl - \W2 • (59) 



This was then used to show that the quantity J (cf. (46)) satisfies the Fisher information 
inequality 

J(Xi Ha X2) < AJ(Xi) + (1 - A) J(X2) . (60) 



The proof of (61 ) follows immediately from the monotonicity (45 ) of the divergence-based Fisher 



information, its additivity (44), the reparametrization identity (43), as well as the compatibility 



properties (58), (59). We will omit the corresponding argument here; it was discovered in the 



classical context by Zamir [36] . 
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Here we argue briefly that the quantity (p2j) satisfies an analogous inequality, that is, 

J(Xi Ha X2\EiE2) < \J{Xi\Ei) + (1 - \)J{X2\E2) . (61) 

Indeed, it is clear that the identities (58 ) and (59 ) still hold if we replace the maps e*^^, £^ ^ ^^ 



and W^ by their 'stabilized' versions (obtained by adjoining an identity) 



^tc 


l-> 


e*^ ® Ie,E2 




l-> 


gX,X,^Y^j^^^^ 


W? 


t-)- 


y^l ® Ie,e, . 



Furthermore, the quantity J{A\B)p^g is also motononous (cf. (53)). Because each of the 



terms J{p ' '']9)\q=q constituting J{A\B)p^g is additive and satisfies the reparametrization 
identities, we can apply Zamir's proof again (carrying along E1E2) and obtain the conditional 



Fisher information inequality (61). 



8.3 The conditional entropy power inequality for Gaussian states 

The entropy power inequality we prove relates the conditional entropy S{Y\EiE2)a of the state 

CTySiEa = (^A ® lEiE2)iPXiEi ® PX2E2) '■= P{XimxX2)EiE2 

to the conditional entropies S{Xj\Ej) of the two (Gaussian) states px e , j = 1,2. 

Theorem 8.1 (Conditional entropy power inequality for Gaussian states). Let pxiEi and PX2E2 
be arbitrary Gaussian states. Then 

S{Xi Ha X2IE1E2) > \S{Xi\Ei) + (1 - \)S{X2\E2) . 

Note that the proof outlined here combined with the discussion in [20] (respectively Stam's 
proof [29]) should also provide the inequality 

S{Xim^X2\EiE2) y ^ S{Xi\Ei)/n , ^ S{X2\E2)/n 

- 2 2 

where Xi and X2 have n modes. We do not discuss this version here for brevity. 



Proof. Let the covariance matrices and displacement vectors of px^Ej , j = 1, 2, be given by (56 ). 
The corresponding covariance matrices are 

Mx,it)E, = (f(f ^^f ) where Mx.^t) = Mx, + tl . 



For t > 0, define the function 



b{t) := 5(Xi(t) m^X2{t)\E^E2) - A5(Xi(t)|Ei) - (1 - \)S{X2{t)\E2) , 
where the entropies are evaluated on the states px{t)E-, j = 1, 2 and the result p(XY(t)my,X2(t))ExE2 



of letting these interact with the beamsphtter (as discussed in Section 8.1). According to the 
'stabilized' version of ((58|, we have 



b{t) = Si{X,mX2m\E^E2) - A5(Xi(t)|Ei) - (1 - A)5(X2(t)|i?2) 
= SiYit)\E,E2) - XSiX,it)\E^) - (1 - X)SiX2it)\E2) . 
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This shows that 6{t) is the difference of conditional entropies of time-evolved states for different 



initial states py E1S2 5 PxiEi and PX2E2 at t = 0. We conclude with Lemma 6.2 that 



lim S{t) = . (62) 

t— ^00 



On the other hand, we have according to the de Bruijin identity (Theorem 7.3) 



4(5' (t) = J(Xi(t) Ha X2{t)\EiE2) - AJ(Xi(t)|^i) - (1 - X)J{X2{t)\E2) 



This identity, together with Fisher information inequality (61 ) imply that 6'{t) < for all t > 0. 



With (62), this shows that 5(0) > 0, which is the claim. D 
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